Approaching the Dirac point in high mobility multilayer epitaxial graphene 
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Multi-layer epitaxial graphene (MEG) is investigated using far infrared (FIR) transmission exper- 
iments in the different limits of low magnetic fields and high temperatures. The cyclotron-resonance 
like absorption is observed at low temperature in magnetic fields below 50 mT, allowing thus to 
probe the nearest vicinity of the Dirac point and to estimate the conductivity in nearly undoped 
graphene. The carrier mobility is found to exceed 250,000 cm 2 /(V.s). In the limit of high tem- 
peratures, the well-defined Landau level (LL) quantization is observed up to room temperature at 
magnetic fields below 1 T, a phenomenon unique in solid state systems. A negligible increase in the 
width of the cyclotron resonance lines with increasing temperature indicates that no important scat- 
tering mechanism is thermally activated, supporting recent expectations of high room-temperature 
mobilities in graphene. 



The quality of electronic systems, measured in terms 
of the carrier mobility, and determined by carrier scat- 
tering processes, is a major issue in material science both 
from a fundamental physics and an applications point of 
view. This issue is now critical for further progress in the 
research on 2D allotropes of carbon 0, 0, [E @1 • In 
exfoliated graphene, charged impurities are recognized 
as the dominant source of scattering in transport ex- 
periments @, S, @| • The presence of charged impurities 
even prevents experiments in the vicinity of the Dirac 
point, due to the electron and hole puddles which in- 
evitably emerge in these samples when lowering den- 
sity [lOj. Recently, the first experiments on suspended 
graphene [111, EH , showed a promising way to come closer 
to Dirac point, but nevertheless, the exact role of indi- 
vidual scattering mechanisms, especially close to Dirac 
point, together with the nature of the minimum conduc- 
tivity remains unclear [H Q E [IE E3 • 

In a parallel track to the exfoliated-graphene research, 
the equivalent problems arise in the study of epitaxial 
graphene. Multi-layer epitaxial graphene (MEG) can be 
relatively easily thermally decomposed from SiC on the 
C-terminated surface and consist of electrically decou- 
pled graphene layers 18, 1^. It is highly conducting in 
the very close vicinity of the interface with SiC, as a result 
of charge transfer from SiC [|, [2pj . The subsequent lay- 
ers are practically neutral 21]. Electrical measurements 
on epitaxial graphene samples probe mostly the low resis- 
tance interface layers. However, the quasi-neutral part of 
the MEG can be investigated by FIR optical absorption 
technique in an applied magnetic field [21 1. 

Both epitaxial and exfoliated graphene have been ex- 
tensively investigated at elevated temperatures. By anal- 
ogy with room temperature ballistic transport in car- 
bon nanotubes [22I ] , very high mobilities are expected for 



graphene at room temperature. This is indeed suggested 
by the investigation of the temperature dependence of 
the extrinsic and intrinsic scattering mechanisms in ex- 
foliated graphene [23| and in the weak temperature de- 
pendence of the resistance in epitaxial graphene [ij]. 

In this Letter, we focus on carrier scattering processes 
in the majority graphene layers with a low carrier concen- 
tration and show their extraordinary quality, expressed 
by the temperature independent carrier mobility in ex- 
cess of 250,000 cm 2 /(V.s). This information is extracted 
from the analysis of inter-LL transitions which are inves- 
tigated in the relevant limits of high temperatures and 
low magnetic fields. 

The investigated graphene sample was prepared by 
thermal decomposition from a 4H-SiC substrate [H and 
contains around ^100 graphene layers. The sample was 
characterized using micro-Raman experiment, which re- 
vealed similarly to Ref. [24| the presence of decoupled 



graphene layers as well as of additional graphite resid- 
uals. The Raman experiment with the laser spot of di- 
ameter ~3 [wa shows an extraordinarily narrow 2D band 
(~20 cm" 1 ) compared to values reported on exfoliated 
graphene 25j . To measure the FIR transmittance of the 
sample, the macroscopic area of the sample (~4 mm 2 ) 
was exposed to the radiation of a globar or mercury lamp, 
which was analyzed by a Fourier transform spectrometer 
and delivered to the sample via light-pipe optics. In ex- 
periments performed at T — 2.0 K, the light was detected 
by a composite bolometer placed directly below the sam- 
ple, while at higher temperatures an external bolometer 
was used. All presented transmission spectra are normal- 
ized by the sample transmission at B = 0. 

The low temperature FIR transmission spectra of the 
investigated sample, shown in Fig. [IJ are typical of the 
optical response of MEG 2l| . Considering the LL spec- 
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FIG. 1: (color online) Part (a): The FIR transmission T 
plotted as — In T as a function of the magnetic field at 
T = 2.0 K. The dashed lines denote the expected transitions 
for c = 1.02 x 10 6 m.s -1 . The inset shows the transmission 
spectrum at B = 0.3 T. Part (b): FIR transmission taken 
at T = 2 K in low magnetic fields. For clarity, successive 
spectra are shifted vertically by 0.05. The part (c) shows the 
peak position, width and area for the L_i(o) — > I J o(i) transi- 
tion. The dashed line in part (c) is a least squares fit to the 
peak positions. 



trum in graphene: E n = sign(n)c^2eHB\n\, four absorp- 
tion lines, denoted by arrows in the inset of Fig.[TJa) and 
showing a VB-scaled blueshift, can be identified as inter- 

>L m with 
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LL transitions L_ m — >L m+ i and L_ 
m = 0,1,2 and 3. The corresponding Fermi velocity 
was evaluated as c = (1.02 ± 0.01) x 10 6 m.s^ 1 . No devi- 
ations from the single particle model due to many-body 
effects are found in MEG in contrast to recent experi- 
ments performed on exfoliated samples 2^, 27 1 . Notably, 
the spectra of this ~100 layer sample present no transi- 
tions symmetric around Dirac point (L_ m — >L m ) which 
are characteristic of bulk graphite [28| . 

Henceforward, we focus on the main line in the spectra, 
Lo(-i) "~ ►Li(o)) which necessarily corresponds to transi- 
tions from or to the vicinity of the Fermi level. Following 
this transition with the magnetic field, we find no devi- 
ation from the \/S-scahng and the line is still visible at 
B ~ 40 mT when it is centered at an energy of «7 meV. 
The FIR experiment, thus allows to probe the very close 
vicinity of Dirac point, hardly accessible in the current 



transport experiments [i , 11 , 12| , and also shows that the 
linearity of the density of states is preserved down to the 
distances of few meV from Dirac point. The disappear- 
ance of the line below B ~ 40 mT gives an estimation 
of the LL filling factor v ss 6 at this field with the cor- 
responding carrier density of no ~ 5 x 10 9 cm~ 2 . The 
possible remanent field of the solenoid of < 5 mT limits 
the accuracy of this estimation to 10%. Note that this 
density is about three orders of magnitude smaller than 
the carrier density measured in transport for equivalent 
samples 0, 0, [29} , but electrical conductance is governed 
by highly doped graphene layer(s) close to the SiC sub- 
strate. While the FIR experiment probes all layers in the 
sample, the strongly doped layer(s) give no contribution 
to the FIR spectra in the presented region of energies and 
magnetic fields. 

The shape of the main line can be well reproduced 
by a simple Lorentzian curve. The results of the fitting 
procedure are shown in Fig. [He), where the peak posi- 
tion, area and width are plotted as a function of \[B. 
Both the peak position and the area show a linear in- 
crease with %/£?, in agreement with expectations for a 
single graphene layer [21( . The latter also suggests a rel- 
atively good homogeneity of the carrier density, as the 
significant presence of more doped regions would result 
in a superlinear rise. An interesting evolution with mag- 
netic field is seen for the linewidth 5E. Starting from 
SE ~ 4 meV at the lowest magnetic field B = 60 mT, 
where the line-shape analysis is possible, the width in- 
creases nearly linearly with \/B up to almost 8 meV at 
1 T. At higher magnetic fields, the increase continues but 
as a sublinear function of \f~B. The observed broadening 
cannot result from the electron-hole asymmetry, where a 
broadening which varies linearly with B is expected [30j j . 
The \/i?-dependent broadening of LLs was suggested by 
Shon and Ando [l3j for the case of both short- and long- 
range scatterers, whose strength is independent of car- 
rier density. Whereas the short-range scatterers should 
induce an identical broadening of all LLs, the n = LL 
should be broadened by a factor of V2 more than other 
LLs for long-range scatterers. As the width of the main 
line is not enhanced compared to the other transitions, 
we can conclude that short-range scattering is probably 
dominant here. 

To compare our results with the recent transport ex- 
periments, a simple estimation of the scattering time, 
mobility and conductivity are presented. From the width 
of the main transition 5E, the carrier scattering time can 
be estimated using r = 2K/5E, also used to compare FIR 
transmission and transport experiments performed on ex- 
foliated graphene [2^,[3l[. The scattering time, obtained 
for our macroscopic sample, is r ~ 300 fs at low mag- 
netic fields (SE s» 4 meV) and decreases with increasing 
B to t ~ 150 fs at B = 1 T (SE » 8 meV). Equivalent 
scattering times r ~ 200 fs are nowadays reported for 
single-flake graphene on a Si/Si02 substrate at densities 
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FIG. 2: (color online) Temperature dependence of the trans- 
mission spectra taken at magnetic fields B — 0.8 and 1.0 T 
in parts (a) and (b), respectively, showing the L-i(o) - > Lo(i) 
transition up to the room temperature. The blue curves rep- 
resent the Lorentzian fits whose parameters are plotted in 
Fig. [3] For clarity, successive spectra in both parts are shifted 
vertically by 0.15. 



around 10 12 cm 2 but this time decreases rapidly with 
decreasing carrier density 
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The scattering time of 
t ~ 260 fs was also reported in charged grap hene layer 
(4x 10 12 cm -2 ) at SiC/graphene interface [29(. Recently, 
scattering times r ~ 100 fs were achieved at density down 
to ~ 10 10 cm -2 in suspended graphene flU H3|. 

The lowest field of B = 40 mT for which the well- 
defined absorption line is observed at huj c w 7 meV al- 
lows an independent estimation of the lower bound for 
t, and the mobility /1. Using the semi-classical con- 
dition lo c t > 1, allowing the carriers to complete one 
cyclotron orbit without scattering, gives a minimum re- 
laxation time t > ft/7 meV~ 100 fs independent of 
but in a good agreement with the previous estimation 
based on the linewidth SE. The condition lu c t > 1 can 
be rewritten as fiB > 1, which leads to the mobility 
fx > 0.25 x 10 6 cm 2 /(V.s). Note that this value repre- 
sents only the lowest bound, as the main absorption line 
disappears from the spectrum due to the complete filling 
(depopulation) of n = 1 [n — — 1) LL and not due to the 
LL broadening. It remains an open question how this mo- 
bility would evolve with increasing carrier density, since 
the scattering mechanisms are different from those in ex- 
foliated graphene, where ionized impurities dominate and 
the mobility remains nearly constant 11, 13, 3, 17|. 

To estimate the conductivity, we can use the standard 
relation a(e) = (e 2 /2)D(e)c 2 t, derived in the Boltzmann 
transport theory. This formula is applicable for graphene 
systems with a homogeneous density of carriers 0, 
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FIG. 3: (color online) Parameters of the Lorentzian fits to the 
transmission spectra in Fig. [2] Red and blue color represent 
data for B = 0.8 and 1.0 T, respectively. Whereas position 
and width of absorption lines in parts (a) and (b) show none 
or a very a weak temperature dependence, the peak area in 
the part (c) significantly decreases with T, nearly following 
the theoretical expectation 1 — 2/i depicted by solid lines. 
The inset schematically shows the LL occupation and possi- 
ble dipole-allowed absorption processes between LLs at finite 
temperature in undoped graphene (Ef = 0). 



fulfilled in our case, but distinctive deviations from a 
presumably more precise self-consistent Born approxi- 
mation are predicted in the immediate vicinity of Dirac 
point 13]. To calculate the density of states D(e), we 
recall that no deviation from the \/B scaling of the posi- 
tion of the main line is observed, and thus the expression 
D(e) — g s g v \e\/ (2irc 2 h 2 ) remains valid down to a few 
meV from Dirac point, where g s = g v = 2 denote the 
spin and valley degeneracy, respectively. Finally, we ob- 
tain the conductivity in the form a(e) = (2e 2 /ft)|e|(r/ft). 
In the limit of a very low magnetic field, when the LL 
quantization energy is comparable to the LL broadening, 
we find SE = 2ft/r — >?« 3 meV and the Fermi level |e| — > 
chy/nno ~ 8 meV which implies a zero-field conductiv- 
ity cr = (Ae 2 /h)(\e\/5E) » 10e 2 /ft, slightly larger than 
the minimum conductivity reported in Refs. [i[ 0, EH, 
but corresponding to conclusions for relatively clean sam- 
ples HEEa. 

Having discussed the spectra taken at T = 2.0 K, we 
turn our attention to measurements at elevated temper- 
atures in Fig. [2j In these data, we can follow the evo- 
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transitions at B — 0.8 and 



1.0 T in the interval of temperatures T = 15 — 275 K. 
Apart from an apparent decrease in the absorption line 
intensity, no other effects are noticeable. Hence, the well- 
defined LLs are uniquely observed in graphene at room 
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temperature and for magnetic fields below 1 T, much 
lower than reported up to now 32]. The blue curves 
in Fig. [3 are Lorentzian fits to the data, and the ob- 
tained parameters, the peak position, width and area are 
plotted in Fig. [3] Based on these results, we can draw 
the following conclusions. The constant position of the 
absorption line shows that the Fermi velocity is indepen- 
dent of temperature. The extremely weak increase of 
the linewidth shows that no important scattering mecha- 
nism is thermally-activated, indicating that the mobility 
|U > 0.25 x 10 6 cm 2 /(V.s) remains practically unchanged 
up to room temperature. For comparison, the mobility 
in standard 2D GaAs systems at room temperature does 
not exceed 10 4 cm 2 /(V.s), being intrinsically limited by 
electron scattering on optical phonons 33] . Scattering on 
neutral centers, independent of carrier momentum and 
in consequence raising the temperature independent mo- 
bility is likely the predominant scattering mechanisms in 
the investigated graphene layers which are well separated 
from the substrate and screened by highly conducting in- 
terface graphene sheets. Such mechanism is consistent 
with the above discussed observation of the nearly \J~B- 
dependence of the linewidth. For comparison, the scat- 
tering on ionized impurities determines carrier mobility 
in exfoliated graphene, which in contrast may lead to a 
temperature dependent mobility at high temperature, in 
case of a non degenerate gas. 

The decrease in the strength of the transition can 
be understood by considering the influence of the ther- 
mal population and depopulation of LL's on the prob- 
ability of the absorption process. The probabilities for 
the transitions Lo — >Li and L_i — >Lo are expressed as 
P ^i oc (f - /i) and P_i^ oc - / ), respec- 

tively, yielding the total probability P cx (f-i — /i), 
where /„ is the occupation of the n-th LL. Assuming 
undoped graphene (Ep = 0), the relative decrease in the 
intensity of the main absorption line can be written as 
7(T) = 1 — 2fi which reproduces the data in Fig. [3fc) 
extremely well. To calculate fi , the optically determined 
Fermi velocity c = 1.02 x 10 6 m.s" 1 was used. Hence, the 
only pronounced effect observed in our FIR spectra with 
increasing temperature is a population effect due to the 
thermal excitation of carriers, as schematically shown in 
the inset of Fig. EIc). 

In summary, MEG has been investigated using LL 
spectroscopy in low magnetic fields and at temperatures 
ranging from 2 to 275 K. Well-defined inter-LL absorp- 
tion are observed at low temperatures down to magnetic 
fields of 40 mT. The LL spectroscopy thus probes the 
electronic states in the immediate vicinity of Dirac point, 
which are hardly accessible in experiments on exfoliated 
graphene llj, [l2j, [lj, lljj due to the inevitable presence of 
electron-hole puddles. The conductivity of nearly neutral 
MEG deduced from our optical data is found to be con- 
sistent with recent theory [9|, . The temperature inde- 
pendent width of the spectral lines, which arise from the 



temperature independent LL broadening, indicate very 
weak thermally activated scattering processes. There- 
fore, the carrier mobility in excess of 250,000 cm 2 /(V.s), 
evaluated from low temperature data, remains practically 
constant with increasing temperature, confirming the re- 
cent expectations of record room-temperature mobilities 
in graphene [23| . MEG is an interesting material to probe 
the properties of Dirac particles in the immediate vicinity 
of Dirac point and has a potential as a platform for high 
performance graphene-based electronics. However, real- 
ization of any electric transport experiments on nearly 
neutral MEG will require passivation of the highly con- 
ducting graphene interface layers. Encouraging research 
along these lines is already underway [34| . 
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